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Vortex matter in mesoscopic two-gap superconducting disks: 
influence of Josephson and magnetic coupling 
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The effects of the couphng between two electronic condensates in two-gap mesoscopic super- 
conductors are studied within the Ginzburg-Landau theory using a finite difference technique. In 
apphed magnetic field, we derive the dependency of the size of the vortex on the sample size and 
the strength of the Josephson couphng. In addition, we elaborate on the dependence of the critical 
temperature and field on the parameters of the coupled condensates. We demonstrate further the 
existence and stability of fractional states, for which the two condensates comprise different vorticity. 
Moreover, we also found pronounced asymmetric fractional states and we show their experimentally 
observable magnetic response. Finally we introduce the magnetic coupling between condensates, 
and study in particular the case where one band is type II and the other is type I, i.e. the sample 
is effectively of I.x type. The calculated M{H) loops show a clear signature of the mixed type of 
superconductivity, which we find to be strongly affected by the ratio of the coherence lengths in the 
two condensates. 



PACS numbers: 74.78.Na, 74.25.Dw, 74.25.Qt 



I. INTRODUCTION 

MgB2 is the first superconductor unambiguously 
shown to possess two superconducting gapSfi Since its 
discovery in 2001;^ a lot of research was conducted on 
this specific material as well as on two-band supercon- 
ductors in general. In its class of binary compounds and 
metallic superconductors, MgB2 turns out to have the 
highest critical temperature known today, Tc = 39 K. Its 
bulk critical field is strongly anisotropic: 3.5 T along the 
c-axis of the crystal and 17 T in the a6-planc and can 
reach as high as 43 T in films i^"— 

While the mechanism of its superconductivity is not 
yet entirely understood, it has been experimentally 
proven that MgB2 has two separate superconducting 
gaps. For example, in Refs.— >2^ the separately imaging 
of the TT- or cr-bands was demonstrated. On the theo- 
retical side, one considered two order parameters to de- 
scribe the superconducting properties of MgB2. One of 
the first Ginzburg-Landau (GL) descriptions of multigap 
superconductors was developed by Zhitomirsky and Dao, 
starting from microscopic theory J^ Fitting to experimen- 
tal results, the authors pinpointed the values of several 
GL parameters relevant for MgB2, and derived analytical 
expressions for the critical parameters. In Ref<^ the same 
authors discussed the anisotropy of Hc2 within the GL- 
framework. They considered only the direct exchange 
of Cooper pairs between condensates, i.e. the so called 
'Josephson' coupling. In the work of Askerzade et al. 
a different interaction between the bands was investi- 
gated - the drag eSectr^^ which is described in the GL- 
formalism through the coupling of the gradient terms of 
the two condensates. The apparent agreement with ex- 
periment led these authors to fitting parameters for the 
GL-modelof MgBa. 

MgB2 is generally accepted to be a typc-II supercon- 
ductor. However, in a very clean sample, Moshchalkov 



et al. estimated that one of the bands could be type-I 
and the other type-IlJ^ii^ The resulting system exhibits 
behavior that cannot be attributed to either type, thus 
the name type 1.5 superconductivity seemed credible. In- 
deed, the authors found a strong clustering of vortices, 
a phenomenon which they ascribe to a combination of 
attractive and repulsive vortex- vortex interaction. Actu- 
ally, Ref.— reported a positive surface energy for vortices 
whenever the coherence lengths of the two bands are com- 
parable. In Refi^ the semi-Meissner state was predicted 
theoretically for a two-gap superconductor and for su- 
perconductors which do not belong to either of the two 
classes type-I or type-II. The possibility of vortices car- 
rying non-integer flux was studied in RefsJ^. 

Surprisingly, virtually all studies done on two-gap su- 
perconductors (TGS) up-to-date concern bulk samples. 
It is known however that mesoscopic superconductivity 
bears a number of fascinating phenomena, ranging from 
specific vortex states, to enhancement of critical param- 
eters by quantum tailoring. The only existing example 
of such a study is the work of Ghibotaru et al. on meso- 
scopic disksfi^ As a novelty, the authors found that frac- 
tional vortex states (when bands have different vorticity) 
can be realized in a TGS and can even be thermody- 
namically stable. However, those results turned out to 
be specific to the case of very weak coupling and not 
realistic for MgB2. 

In this paper we analyze the fundamental properties 
and vortex matter of mesoscopic disk-shaped two-gap 
superconductors using the Ginzburg-Landau formalism, 
where the electronic exchange between condensates oc- 
curs through Josephson coupling, and magnetic exchange 
between condensates is allowed for. The latter mecha- 
nism has not yet been studied in detail up to now. The 
paper is organized as follows. In Section II, after describ- 
ing the theoretical approach, we focus on the effects of 
Josephson coupling on the the size of a vortex core, the 



unique vortex states and their H—T stability regions, and 
the critical temperature and field as a function of cou- 
pling strength. In Section III, we introduce the screening 
of the magnetic field into the theoretical formalism, and 
illustrate the influence of the magnetic coupling between 
the two condensates on the vortex states, particularly in 
the case of type I.x superconductivity. Magnetic signa- 
tures of the different features are discussed in the light 
of potential observation by magnetometry Finally, our 
findings are summarized in Sec. IV. 



II. JOSEPHSON COUPLING 

A. Theoretical formalism 

It is widely accepted that the high critical tempera- 
ture of MgB2 arises due to the coupling of the super- 



conducting bands which effectively reinforce each other. 
However, the exact nature of the coupling is not fully 
understood, and possible scenarios are the exchange of 
electrons. Cooper pairs, interaction between the respec- 
tive supercurrents, interaction through the internal mag- 
netic field, etc. Microscopic ab-initio calculations have 
not been able to pinpoint the key interaction. In what 
follows, we will consider the Josephson coupling between 
the bands, resulting from the tunneling of the Cooper 
pairs from one band to another. This is incorporated in 
the Ginzburg-Landau (GL) energy functionali^ through 
an interaction term dependent on the order parameter of 
both bands and proportional to F, the Josephson cou- 
pling strength: 



J 



AF 



E 



2m„ 



{^ihV - —A)* 



aJ'^r. 



:/3n|^. 



r(*t*2 + *2*i)) 



dV, 



(1) 



where a„ ~ Q:no(l ^ T /Ten) and /3„ are the GL coeffi- 
cients, and n is the band index. This results in a set of 
nine parameters describing a two-gap system: aio, 0:20 j 
/3i, /32, the Cooper-pair mass mi and m2, F, and the 
critical temperatures Td and Tc2- While each band has 
its own intrinsic critical temperature, Josephson coupling 
causes both bands to survive up to a higher critical tem- 
perature, Tc > max(Tci,Tc2). While the GL functional 
is derived for T < Tc, experience with mesoscopic single- 
gap superconductors indicates that the GL equations in 
practice are valid much deeper into the superconducting 
state. In Eq. ([T]) we then introduce temperature indepen- 
dent units, in order to rewrite it in a dimensionless form. 
We express the free energy of the system in units of _Fio = 
"lo/Z^i^ length in units of ^10 (^„o = fi/V^ 
vector potential in Aq ~ hc/2eS,io, 
ters in *„o = ^«o(r = 0,F = 0,if = 0) 
and the temperature in Td . This reduces the set of nec- 
essary parameters to five: S = ^10/^20; " — ^,10/1,20' 
m = mi/m2, T^r = Tc2/Tci and 7 — T/aio- In addition, 
we have two external tuneable parameters: the tempera- 
ture T and the applied field H. The flux is defined as 
the externally applied flux. We flrst consider an extreme 
type-II case, and neglect the self-induced magnetic field 
in the sample. 



The minimization of the energy functional leads to 
the two-band GL equations. After the scaling described 
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above, the equations for the order parameters read: 

(-zV - A) Vi - (1 - T - IV'i nVi - JV'2 = 0, 

-(-ZV - Afi'2 - f 1 - 7^ - IV'2p) ^-2 - — V'l = 0. 
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This system of non-linear coupled differential equations 
we solve numerically on a square grid of typically 128x128 
points. The details of this procedure can be found in 
Ref^. 

In the following analysis we neglect the screening of 
the magnetic field. This is justified for an extreme type- 
II material, or any sufficiently thin sample. An applied 
vector potential A = {^Hy,—^Hx,0) results in a mag- 
netic response resulting from a total supercurrent: 
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The free energy functional in dimensionless units reads: 
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Let us here address several direct implications of 
Josephson coupling. It is clear from Eq. ([T]) that the sign 
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FIG. 1: (Color online) (a) The dependence of the observed 
vortex size Rv on the number of grid points in the numerical 
mesh N . (b) The relative error in found vortex size vs. grid 
spacing h, for two different sizes of the superconducting disk. 



of 7 determines the relative phase shift between the or- 
der parameters in the two condensates - either « when 
7 > or « TT when 7 < - in order for the coupling 
term to provide a negative energy contribution. How- 
ever, the sign of 7 has no influence on observables such 
as the Cooper-pair density and magnetic response of the 
sample. The general consequence of 7 coupling is an in- 
jection of Cooper pairs from one band into the other and 
vice versa, thus increasing the stability of the supercon- 
ducting state. In other words, the average Cooper-pair 
density always increases with 7. In the absence of an 
applied field, the ratio x ~ "01 /'/'2 can be found from 
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analytically derived from the GL-equations. From this 
we find that in the limit 7 — >■ c» a constant ratio 
V'i/V'2(H=0) = yJ^/ma/5 is reached, independent of 
temperature. 

The next section is dedicated to an analysis of the size 
of the vortex core. Before we get into the physics of 
the problem, we here address some numerical issues fol- 
lowing from mapping of the superconducting disk on a 
square numerical grid. For obvious reasons, the influ- 
ence of the resolution of the numerical grid A^, on the 
observed vortex size Ry, is significant. In Fig. [TJa) we 
show the Rv vs. N (definition of Rv is given in the 
next section). With increasing grid density, the numeri- 
cal error decreases, and the vortex size converges towards 
the i?o value with dependence Rv = -Ro + h/N . By a 
fitting procedure i?o and b can be determined and the 
corresponding curves are represented by the solid lines 
in the figure. The dashed lines are the asymptotes with 
value i?o- In Fig. [TJb) the relative error is plotted as a 
function of ft., the grid spacing (proportional to sample 
size and inversely proportional to N) . We conclude that 
due to numerics the vortex size is always slightly over- 



estimated, with overshoot increasing for smaller samples 
(where influence of the boundaries is more pronounced). 
Nevertheless, with a resolution of 10 points per coher- 
ence length we get a relative error under 5%. Although 
higher grid density obviously improves the results, we re- 
frain from using a density above 10 points/^, in order to 
optimize the speed of the calculation. 



B. Size of the vortex core 

The coherence length is the characteristic length scale 
over which the order parameter changes. It is there- 
fore intuitive that the size of the vortex core is propor- 
tional to the coherence length in bulk superconductors. 
In single-gap materials, the coherence length is propor- 
tional to l/^/l — T/Tc in the temperature range where 
the Ginzburg-Landau theory is valid. Here we show 
that in two-gap superconductors the coherence length is 
strongly affected by the coupling parameter 7, generally 
in an opposite manner from temperature. For compari- 
son, we can use the coherence length obtained from the 
expression for the second critical field derived in Ref.— , 
through the relation Hc2 = *&o/27r^^, with $0 being the 
flux quantum. Deviations are a priori expected, since al- 
ready experiments of Refs.— >2^ found a discrepancy be- 
tween the vortex size and the coherence length deduced 
from the second critical field. 

To estimate the coherence length, wc will numerically 
determine the size of the vortex core, a quantity which 
is not uniquely defined. In the following calculations, we 
examine the single- vortex state in a MgB2 superconduct- 
ing disk exposed to a field providing three flux quanta 
through the sample. For the definition of the vortex size, 
two possibilities are considered in the literature: (i) The 
vortex size is determined by the distance from the center 
of the vortex to the contour where the Cooper-pair den- 
sity (CPD) recovers to some percentage of its maximal 
value in the sample, denoted as Rv,cpd', and (ii) The 
vortex size is the distance from the center of the vortex 
to the first contour where the supercurrent js reaches its 
maximum, denoted as Ry.j ■ The problem of the first def- 
inition is the arbitrary threshold value for the criterion, 
but also the fact that we have two Cooper pair densi- 
ties, which makes the single vortex size ambiguous. As a 
threshold we take 80 %, since this allows for a more pre- 
cise vortex size determination, and we will consider only 
the first condensate. On the other hand, the second defi- 
nition involves coupled condensates and thus provides us 
with a unique vortex size. We therefore adopt the second 
definition to describe the vortex size in the rest of this 
work. Contrary to the bulk case, in our mesoscopic disks 
both definitions render a vortex size dependent on the 
radius of the disk Rjj, as vortex currents in the center 
of the sample can interact with Meissner currents decay- 
ing from the edge inwards. However, while Rv,cpd in 
each band saturates for Rjj — > od, this is not the case 
for Rvj'- it develops a linear dependence on Rd. We 



extracted the exact dependence of both definitions of the 
vortex size Ry on the disk size Rd, which led us to a 
universal formula (valid for both definitions, but with 
different coefficients) : 
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where Ryo is the vortex size independent of the sam- 
ple size, c is a length coefficient and h is the slope of 
Rv vs. Rd for large Rd- From fitting of our numerical 
data, collected at different T, ^, a and m, we obtained 
{c = 1.90, /i = 0} for Rv,cPD and {c = 3.40, h = 0.006} 
for Rv^j ■ These coefficients are valid for single gap super- 
conductors as well as for two-gap superconductors, even 
when the coherence lengths of the two condensates are 
very different (e.g. for small a, see Fig. [2]). The function 
gives an excellent estimate of the vortex size for to 7^ 1, 
for disks larger than 5^. As shown in Fig. [2l deviation 
from the given function does occurs for specific choices 
of m and a, especially for small disks, but the relative 
error remains under 5%. 

Both previously given definitions of Ry are illustrated 
in Fig. [3l When disks are too small, the vortex size is 
not always fixed by the disk size only, e.g. when T > Tcr 
and 7 < 0.1, when the coherence lengths differ much and 
the interaction of the vortex with the Meissner currents 
becomes too different in the two bands. We found how- 
ever that when Rd > lO^io these mesoscopic effects have 
only a minor influence and the correspondence between 
Rv and Rd becomes predictable again. 

With the established dependence of the vortex size on 
the size of the sample, we can more precisely determine 
the actual influence of the coupling 7 on the vortex size. 
In particular we will look at the behavior of Rvo.j as a 
function of T, Tcr and 7 (parameters a and m remain 
fixed at realistic values for MgB2). 



°°°^.^^^'°'°""^rVh 



So5HB0tiflfla 



.82. 



_■■■ oo„[ 

, ■ o a a S "„ Q^ 



□ u g B B °-ai 



» 



.-...-.8^0 
, °oo8°° , 

JooO° n' 

o 00 o " ■ 



5d 



■ 


a=0.2 


m=5 




a=0.2 


m=15 




a=0.2 


m=1 


n 


a=0.85 


m=1 



10 



15 



20 



25 



30 



FIG. 2: (Color online) The relative error of the result of Eq. 
([5]) compared to the observed vortex size, as a function of 
disk size, for indicated different values of parameters m and 
a, each with three different combinations of 7 and T in order 
to cover as large as possible parameter space in the analysis. 



In Fig. |4]thc size- independent vortex size Ryo is plot- 
ted versus 7, for different temperatures T and Tcr- The 
dots represent the result of our simulations. The gen- 
eral behavior can be described by the following observed 
trends: i) Increasing T induces an increase of the vortex 
size whereas increasing 7 has the opposite effect, ii) De- 
viation from the latter monotonic behavior occurs when 
T sa Tcr and the coupling is weak, see e.g. the curve 
at T = Tcr = 0.4. The reason for the initial positive 
slope is that the second band is revived by the presence 
of coupling but retains its own character (i.e. a larger 
coherence length) since coupling is still weak, iii) For 
T > Tcr, curves with different Tcr but identical T merge 
at 7 = since then only the first condensate survives 
and fully determines the vortex size. In Fig. IH the solid 
curves represent an estimate of the vortex size based on 
the general relation between the coherence length and the 
upper critical field in a single gap bulk superconductor 

Hc2 cx 1/e, 
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g+ia,T,Tc, 



.g_(a,T,r,,)2+4: 



(6) 



with g±{a, T, Tcr) = 1— T±a(l — 7^ ) , based on the 



analytical expression for the critical field of a bulk two- 
gap superconductor. 




taken from Ref.— . We find that the vortex size in our 
samples scales with the coherence length as Ryo = 1.78^, 
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FIG. 3: (Color online) Vortex size, determined through the 
decay of the Cooper-pair density Rv-cpd (a) and determined 
through the maximum of encircling currents Ry-j (b), as a 
function of the radius of the superconducting disk Rd - The 
solid dots represent the numerical results, while solid lines 
show the fit using Eq. O The dashed lines indicate Rvo, the 
fitting parameter corresponding to the sample-independent 
vortex size. 
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FIG. 4: (Color online) Rvo as a function of 7 and tempera- 
ture, for 7 < 1 and for different values of Tc2- a = m — 1. 



which we use to plot the curves in Fig. U) These theoreti- 
cal curves coincide rather well with the data for 7 > 0.25. 
The reason for this is that, when coupling becomes suf- 
ficiently strong, both order parameters tend to have a 
similar spatial distribution and thus also exhibit a similar 
vortex size and coherence length. A good correspondence 
between the data and the fitted curves is also found for 
T ^ Tcr, i.e. when the second condensate exists solely 
due to the coupling to the first condensate, or in the case 
of weak coupling and the second condensate is almost 
depleted, so that it does not influence the vortex size. 
Two regions of discrepancy include T « Tcr {Rvo be- 
haves non-monotonic), and T <^ Tcr and weak coupling. 
For the latter case, the formula still predicts Rvoil = 0) 
to be independent of Tcr, while this is clearly not the 
case. In this regime the vortex size is found to behave 
more like that of a single gap superconductor, but with a 
different critical temperature. By fitting we determined 
a function that describes the behavior of the vortex size 
accurately in this regime as 
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This equation is generally applicable, and effectively 
shows our initial premise that 7 has an opposite influ- 
ence to T. 

We notice however that Eq. ^ also contains the de- 
pendence on a and m. However, this formula can not 
adequately describe the vortex size for m, a much differ- 
ent from 1. since the properties of the two condensates 
can no longer be described by a single coherence length. 
In general we can state that the relation between the 
coherence length and the critical flcld does not hold any- 
more when the individual coherence lengths differ too 
strongly. In Fig. [5] we plot the numerically obtained Ry 
as a function of m, for a small parameter a (with thus 
an acute difference between the coherence lengths in the 
two condensates). The analytic estimate of Eq. (jH]) is 
monotonously increasing with m in this case, and is ob- 
viously not useful for comparison with non-monotonically 



evolving curves in Fig. [5l 

Let us first analyze the limiting case of extremely large 
m. Following from Eq. pi A[) . the second condensate de- 
couples from the first one in this limit. At the same time, 
as seen in Eq. ^, the influence of the second condensate 
on the total current in the system increases. Therefore, 
the size of the vortex Rvj will be fully determined by 
the second condensate (and its nominal coherence length 
^2(2^)), provided that the temperature is below Tcr- Oth- 
erwise, the vortex size is determined solely by the flrst 
condensate (and ^i(T)), since the coupling between the 
condensates is entirely suppressed and the second con- 
densate fully depletes. 

This helps us to understand the behavior of the vortex 
size as a function of m, shown in Fig. [5] At T < Tcr, the 
initial increase of m decreases the coupling of the second 
condensate to the first (to has an opposite effect from 7, 
see Eq. (|II Ap . This causes an increase of the apparent 
vortex size, due to much larger coherence length of the 
second condensate (a = 0.1), and at the same time the 
increase of m makes the supercurrent of the second band 
stronger and thus more deterministic for the magnetically 
detectable size of the vortex (see Eq. ^. Above Tcr 
these two effects become competing, since the first will 
deplete the second condensate and therefore reduce its 
influence while the second enhances the influence of the 
second condensate. These competing effects result in the 
non-monotonic behavior of the vortex size vs. m in Fig. 
[5] At low 771, the influence of the large coherence length in 
the second condensate dominates, whereas at large m the 
coupling disappears and the second condensate depletes. 
At very large 777, all curves for T > Tcr saturate to the 
same value, namely the size of the vortex core in the first 
condensate, in the absence of a second one. 

In Fig. [S] we demonstrate some peculiarities of the 
dependence of the vortex size on the parameter a. a 
was swept down from 1 to 0.005 in a disk with param- 
eters 777 = (5 = 1, Tcr = 0.5 and R/o = 4^io- In the 
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FIG. 5: (Color online) The apparent vortex size in a disk 
with radius K = 30^ as a function of vn, for different 7 and 
T . The ratio of coherence lengths in two condensates is fixed 
at ^20/Cio = 3.162, i.e. a = 0.1. 



absence of coupling, this sweep increases ^20 while ^10 is 
kept constant, as the length unit of the GL equations. In 
the presence of coupling, both the resulting ^2 and the 
resulting ^1 will be influenced (i.e. the Cooper-pair cor- 
relation length in each of the condensates, different from 
the nominal coherence lengths in each condensate sepa- 
rately). Intuitively, one expects that coupling causes vor- 
tex cores in the two condensates to have similar behavior, 
and tend towards similar sizes; instead, for decreasing a 
at temperatures T = 0.45 and T = 0.7 an increase of Ry2 
is observed while Ryi decreases! For stronger coupling 
this effect becomes even more prominent, compared to 
the vortex size at a = 1. At lower temperatures (shown 
for T = in Fig. ©, the behavior of Rv,M) for T = is 
in better concordance with the intuition: Ryi increases 
as Rv2 increases, with the effect growing with coupling. 
However this effect is weak. 

To better understand the behavior of the vortex sizes 
we will invoke the full free energy expression of Eq. ([3]) , 
which can be written d& F = Fi+ F2 + F12, where Fi de- 
pends only on ipi and F12 is the Josephson coupling term. 
The key point here is that a decrease of a will lower F2 
directly, therefore giving more weigth to the other terms 
Fi and F12, i.e. changing a reorders the hierarchy of the 
terms. A stronger F12 stimulates an increase of "01 and 
"02 ■ Since F2, which regulates the size of "02, becomes 
less important, 02 will increase much faster than i/ii. A 
similar behavior of the tp and Ry curves hints to the link 
between the two variables. For an infinite superconduc- 
tor, sweeping a to zero would cause both tpi and 02 to 
diverge. However for a finite (mesoscopic) superconduc- 
tor, the vortex size will eventually exceed the disk size. 
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FIG. 6: (Color online) Top: Ry in each condensate as a func- 
tion of a, for different 7, T. Bottom: Ratio of condensation 
energies of the two condensates (logarithmic scale) vs. a. 
When equal to 1, condensates influence each other with equal 
strength since m = 1. 



thereby effectively suppressing the order parameter and 
preventing the divergence. 

The decrease of Ryi with decreasing a can be ascribed 
to the increase of 0i . However, in the left panel of Fig. [S] 
where temperature is zero, we notice the subtle increase 
of Rvi with decreasing a. Here -02 is still relatively large 
(because of the low temperature) which causes 0i to feel 
a strong influence from 02 due to coupling. As a result, 
for low T, vortex sizes will tend to be similar in the two 
condensates, as we intuitively predicted. In the central 
and the right panel of Fig. [SI i.e. for higher tempera- 
tures, the influence of 02 is reduced, and the increased 
order parameter -01 prevails in determining the vortex 
size - thus Ryi decreases. In the right panel T > T2, 
and the second condensate would be depleted in the ab- 
sence of coupling. This creates yet another regime of the 
Rv (a) dependence: initially Rv2 now decreases with de- 
creasing a from I. The second condensate is completely 
dependent on the first one, and therefore obeys its shape. 
When a further decreases, ^02 is less modified and an in- 
crease of Rv2, is recovered. 
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FIG. 7: (Color online) The size-independent vortex size Rvo 
as a function of temperature (for fixed Td and Tc2), for several 
values of 7. Rd ~ lO^o, m = 1, and a = 1. 

Finally we point out one more interesting artefact. In 
Fig. [7] we show the calculated vortex size as a function 
of temperature, for very weak coupling, i.e. small 7, and 
m = a = 1. The observed kink corresponds to the critical 
temperature of the second band, and smears out when 7 
is increased. The {7 = 0, Tcr = 0.4} curve starts at 
the same value as the {7 = 0,Tcr = 1} curve, since at 
T = there is no dependence of vortex size on Tcr- This 
behavior is observable by magnetic-force, scanning Hall 
probe, or scanning tunneling microscopy, and we expect 
its experimental verification. 

For the {7 = 0,Tcr = 1} curve, the superconducting 
state and supercurrents in both condensates are identical. 
Therefore the vortex size defined on separate condensates 
as well as on the combined system will be equal. For this 



reason the curve coincides with the {7 = 0,Tc2 = 0.4} 
case for T > Tc2 = 0.4, since then the second band is 
depleted and only the first band superconducts. 



C. H-T phase diagrams 

In increasing magnetic field, more vortices penetrate 
the superconducting system. It is known that the sym- 
metry of the vortex states is strongly affected by the sym- 
metry of the mesoscopic sample, as detailed in Refs.— . 
With increasing temperature, the symmetry of the sam- 
ple is even stronger imposed on the vortex matter, and 
it is therefore no surprise that in mesoscopic disks most 
vortex state configurations collapse into a giant- vortex at 
high temperature. We can construct an H — T diagram 
for mesoscopic samples, indicating the area of stability 
of states with different vorticity. Two-gap systems make 
there no exception, but do comprise several particulari- 
ties. In Fig. IS] we displayed the full stability regions of 
all possible vortex states with vorticity i < 7 in a super- 
conducting disk with parameters {Tc2 = 0.44, S = 1.33, 
a = 0.844, m = 1},— which is very similar to MgB2 
except for the coupling parameter, where we took sig- 
nificantly smaller 7 = 0.01. This choice provides more 
complexity to the vortex states, as it allows for different 
vorticities and vortex arrangements in the two bands. 

Indeed, one difference from the single-gap supercon- 
ducting disks is directly visible in Fig. [3 where the 
stability regions of composite vortex states are mush- 
room shaped. In other words, with increasing temper- 
ature, one can exit the stability range of a particular 
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FIG. 9: (Color online) Illustration on how the stability flux- 
temperature regions of integer (a) and fractional states (b) 
arise from the single-gap picture. 



L-state, but then find it again at higher temperatures. 
This shape has the following origin; At high tempera- 
tures (Tc2 <^ T < Tci) the second gap would be com- 
pletely depleted if it wasn't for the coupling. In other 
words, the second band depends completely on the first 
band, and therefore has the same behavior and features 
like the H(T)-boundary. However, for low temperatures 
(T < Tc2) the second band is still active and retains its 
own character, and therefore the stability region bound- 
ary follows quite closely the single-gap stability region. 

In the mushroom-shaped areas, the vorticity in the two 
bands is the same. However, in the shaded areas we found 
vortex states where the vorticity differs from one band 
to the other. As a consequence, the overall, apparent 
vorticity of the sample is no longer integer! These are 
the so-called fractional vortex states. 

As clearly shown in Fig. ^a.), the full H — T stabil- 
ity region of an integer flux vortex state is related to the 
union of the H ~ T stability regions for the given state 
in the two corresponding single-gap condensates. On the 
other hand, the fractional states are found at the inter- 
section of two corresponding single-gap stability regions 
[sec Fig. igb), for the (Li,£2) = (4,5) state]. 

With increasing coupling parameter 7, the vortex 
states in the two condensates are linked together, and 
moreover reinforce each other. In Fig. [10] we show the 
stability region of the integer flux L = 4 state, for three 
values of 7, where the H—T stability region grows with 7. 
We conclude that increasing 7 stabilizes the integer flux 
states, but for the same reason destabilizes the fractional 
states. We discuss the latter further in the following sec- 
tion. 



FIG. 8: (Color online) The magnetic flux-temperature stabil- 
ity regions for different vortex states (vorticity L) in a two- 
band disk of size Rd/£,iq ~ 4, and with Tc2 ~ 0.44, S — 1.33, 
a = 0.844 and 7 — 0.01. In color-coded areas, the vortex 
state is fractional and cannot be represented by a single L, 
but rather as (Li,L2) state, where vorticities in two conden- 
sates are given respectively. 



D. Fractional vortex states 

The existence of fractional states, i.e. states with dif- 
ferent vorticity in the bands, depends strongly on the 
coupling. They survive only at weak Josephson cou- 
pling between the bands, while only integer flux states 




FIG. 10: (Color online) The stability region of the L — 4 inte- 
ger flux vortex state, for three different strengths of coupling 
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FIG. 11: (Color online) The stability regions in (^ — 7 pa- 
rameter space of fractional states with different vorticities in 
the two bands. Parameters of the sample are Rd/ £,10 = 4, 
a = S = m = 1 and T = 0. In insets at the top of the 
figure, we superimposed the logarithmic plots of the Cooper 
pair density in the two gaps on each other (red/blue shades 
for condensates 1/2 respectively) for states indicated in the 
phase diagram by the red numbers. 



are possible at large 7 values. This is illustrated in Fig. 
[TTl where we show that the region of stability of frac- 
tional states shrinks with increasing coupling, but also 
that lower vorticity fractional states are more resilient 
to 7. Another interesting aspect of fractional states is 
their strong affinity to asymmetry. In both condensates 
vortices attempt to form a symmetric shell, but due to 
coupling and different respective number of vortices, the 
final state becomes asymmetric in most cases. For that 
reason, the asymmetry is more apparent at larger cou- 



pling 7. We show several examples through the log-plots 
of the Cooper pair density of the chosen states in Fig. 
rm Note that the fractional state not necessarily con- 
tains vortices in both condensates; for example, inset 1 
in Fig. [11] is the (0, 1) state. Due to coupling, the total 
energy is minimized when regions with depleted order 
parameter in two condensates are on top of each other. 
As a result, the vortex of the second band is attracted to 
the boundary of the sample, where the circulating Meiss- 
ner currents strongly suppress the order parameter in 
the first band. Inset 2 is the (1,2) state for 7 = 0, i.e. 
the condensates are decoupled. This fractional state is 
therefore two-fold symmetric, but when we increase 7 we 
enhance the asymmetric (1, 2) state, as shown in inset 3. 
One vortex of the second band is attracted to the vortex 
of the first band, and the other is attracted to the edge 
of the sample. Finally we show in inset 4 the (4, 5) state, 
at the verge of its stability region, showing maximally 
pronounced asymmetry. Four vortices in both bands sit 
on top of each other, and the remaining, fifth vortex of 
the second band, breaks the symmetry and is gradually 
pulled out of the sample. The found states look similar 
to what was found earlier for Coulomb bound classical 
particles;^ although underlying physics is very different. 

A two-gap mesoscopic system is a prime example of 
a vortex system with competing interactions. Besides 
the vortex-vortex interactions in each band, one must 
take into account the coupling between order parame- 
ters across the bands, and the mesoscopic effect of the 
compression of vortices to the interior by the circulating 
Meissner current that is maximal at the edge. For ex- 
ample, consider the (0, 1) state, where an outward force 
originates from the coupling between the vortex in the 
second condensate and the suppression of superconduc- 
tivity at the edge of the first one. However, this action 
competes with the inward force exerted by the Meissner 
current. This purely mesoscopic effect leads to a tune- 
able position of the vortex in this fractional state: while 
the Meissner current is roughly the same at a given mag- 
netic field, the changed coupling between the condensates 
brings the vortex further to the boundary. This is shown 
in Fig. [1^ as a transition from a fractional (0,1) vortex 
state to an integer L = vortex state with increasing 
coupling. 

As fascinating as they are, the fractional states are 
difficult to find in the ground-state. For example, when 
the coherence lengths of the two bands are the same, 
then the energy landscape in both bands - considered as 
separate single-gap superconductors - are proportional, 
i.e. Fi = a6'^/mF2. All possible vortex states thus have 
their ground state in the same phase space region. The 
total energy of the system, F = Fi -I- i^2, will therefore 
be proportional to the single-gap energy with as a direct 
consequence that fractional states always will have higher 
(or equal) energy compared to the integer states. To 
realize fractional states as the ground state, one therefore 
needs to make the discrepancy between the coherence 
lengths as large as possible. This can be done by taking a 
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FIG. 12: (Color online) Calculated free energy of the frac- 
tional (0, 1) state, as a frmction of Josephson coupling between 
the bands. Insets show contour plots of the Cooper-pair den- 
sity in the second gap, illustrating how the asymmetry grad- 
ually increases with 7 for taken MgB2 parameters and size of 
the disk Rd = 5^io. The 7-sweep was done along the vertical 
dashed line in Fig. 1131 
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FIG. 14: (Color online) (a) The critical field of the mesoscopic 
disk Her as a function of the mass ratio in two condensates. 
Parameters used are Rd ~ 4^io, Tc2 = 0.5, a — 0.1. (b) Idem 
but now as a function of a. 



the case of a fractional state. We will revisit this point 
in the section devoted to magnetic coupling. 



significantly different from one, or by taking temperature 
close to Tc2, when Td > Tc2- In Fig. [13] we show the 
stability and ground state regions of the non-composite 
states in a disk of size R = 5^io at temperature T = 0.4 
(and a = 0.5). 




FIG. 13: (Color online) The stability regions and the ground 
state (colored) regions in (^ — 7-space of the fractional states. 
Taken parameters are Ro/iio = 5, a = 0.844, 5 = 1.33, 
m = 1, Tc2 = 0.44 and T = 0.1, corresponding to MgB2. 



E. The superconducting-normal phase boundary 

As mentioned in preceding sections, in Refs i^^'^^ an ex- 
pression was derived for the second critical field of a bulk 
two gap superconductor, given by Eq. ([7]). That expres- 
sion contains the dependence of the critical field not just 
on 7, but on m and a as well. We have shown that that 
dependence does not describe the vortex properties in 
the two-band samples, but at this point we check its ap- 
plicability for the estimation of the upper critical field of 
mesoscopic two-band disks at a given temperature. Our 
results for the dependence of the upper critical field on 
the mass ratio in two bands are shown in Fig. I14f a). and 
demonstrate perfect agreement with Eq. ([7]), provided 
that the found critical field is scaled by its value at zero 
temperature and in absence of coupling. 




These asymmetric states can be observed in meso- 
scopic two-band samples. As main candidates for such an 
experiment, we select the imaging of only the 7r-band, as 
was done recently in Refi^. Alternatively, scanning Hall 
magnetometry or magnetic force microscopy can both re- 
veal the asymmetric magnetic response of the sample in 



FIG. 15: (Color online) (a) The critical temperature T^r of a 
two-band mesoscopic disk as a function of Josephson coupling 
7, in absence of magnetic field, (b) The upper critical field 
Her vs. 7 for T = 0. Dots represent the numerical data, and 
the solid line is the result of Eqs. ([3[5]). 
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FIG. 16: (Color online) The superconducting-normal phase 
H{T) boundary for the sample with parameters Rd = 4^o, 
5 = 1.33, a = 0.844, Tc2 = 0.44 and m = 1, for different 
values of the Josephson coupling strength. Different colors of 
the dots mean different vorticities. Dashed lines are obtained 
from Eq. with a prefactor described in the text. 

Eq. Q is further applicable for the estimation of 
the H — T supcrconducting-to-nornial phase boundary. 
Namely, equating that expression to zero gives the ex- 
pression for the critical temperature of the two-band sam- 
ple: 



1 



1 + T„ 
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(1 - Tcrf + 4^r, 



ma 



(9) 



This means that the critical temperature of a two gap 
superconductor is always equal or higher than the sum 
of the critical temperatures of the two bands, in the case 
as if there was no coupling. This observation is in contra- 
diction with findings of Ref.— , where it is claimed that 
also a lower Tc is possible, depending on the parameters. 
The authors obtained these results from a microscopic 
derivation. 

Above expressions were originally derived for bulk sam- 
ples. It is already known that the upper critical field 
in mesoscopic superconductors is higher than in bulkj^ 
and it is therefore intuitively clear that Eq. ([7]) would 
not work for the case of two-gap mesoscopic disks. In 
Fig. [15] we show the numerically obtained critical tem- 
perature Tcr and upper critical field Her (corresponding 
to bulk Hc2) versus 7 in disks of size R = 4^io. We 
found that both the dependence of critical temperature 
and field on 7 obey the dependencies given in Eqs. ([Tlini), 
provided that the critical field is scaled to its value in the 
absence of coupling and at zero temperature. In Fig. [!§] 
we show the calculated Tc{H) boundary for different cou- 
pling strengths. Eq. ([7|) can also be inverted to describe 
the dependence of T^ on the applied field. Although de- 
rived for bulk, we find that latter equation nicely fits the 
-ff(r)-curves in Fig. [15] for a mesoscopic disk, after the 
aforementioned scaling of the magnetic field. 



III. MAGNETIC COUPLING 

In the previous section, we assumed the existence of a 
Josephson coupling between two superconducting bands, 
but we neglected the screening of the magnetic field. In 
applied magnetic field, the magnetic response of a two- 
band superconductor follows from the induced supercur- 
rent: 
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Conventionally, the demagnetization and screening ef- 
fects in mesoscopic superconductors are expressed 
through the Ginzburg-Landau parameter k, being equal 
to the ratio of penetration depth A and coherence length 
^. For that reason, we reformulate the equations of sec- 
tion (|II Ap to introduce K2, the GL parameter of the sec- 
ond condensate instead of the parameter (5, the ratio of 
the order parameters in two bands. We start from the 
definitions of ^ and A: 
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which we then substitute in the GL equations to obtain: 
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This form of two-band Ginzburg-Landau equations is 
particularly convenient for comparison with the conven- 
tional types of superconductivity. In the single-gap bulk 
samples, the value of n above or below l/-\/2 determines 
the superconductor being of second or first type, respec- 
tively. For a two-band sample, this distinction is much 
more difficult to establish, since Eqs. (|12m4[) show the di- 
rect infiuence of not only ki and K2, but also the Joseph- 
son coupling and the squared ratio of coherence lengths 
in two condensates a. 

Eq. (|T4l) also shows that two bands are directly coupled 
through the screening currents, and this type of coupling 
we refer to as magnetic coupling. In Fig. [T7] we show 
the calculated magnetization of the disk with radius R = 
10^10 as a function of the applied field (in absence of 
Josephson coupling) , to illustrate how the magnetic field 
couples the two bands. We observe: (i) the magnetization 
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FIG. 18: (Color online) The magnetic response of the frac- 
tional (0,1) and (1,2) vortex states in a mesoscopic two gap 
superconducting disk with parameters Rd = lO^o, ^i = 10, 
K2 = 2, Q = 0.3, 7 = 0.01 and T = 0. The applied field in (a) 
is 0.04Hc2 ((^/(^o = 2) and in (b) 0.0877^2 (<^/<^o = 4). The 
dashed line shows the sample boundary. 



FIG. 17: (Color online) Magnetization versus applied mag- 
netic field {M{H)) loops for a two-gap disk with radius 
Rd = 10^10, and two condensates as single-gap samples, 
in the absence of Josephson coupling. Magnetic coupling 
is included in the calculations, with parameters ki = 3.68, 
K2 = 0.66, and a = 0.06. 



of the coupled system is somewhat higher than the sum 
of the two uncoupled systems; (ii) the found cusps are 
wider, since the Meissner currents in one band screens 
the field in the other band and flux entry is therefore 
more difficult. The net field in the interior of the sample 
is therefore lower, and the critical field is increased. This 
is in accordance with existing experimental findings on 
MgB2. 

However, the information available in literature is also 
often confusing. For example, we identified two sets of 
parameters, both believed to be correct for MgB2. From 
Refs i^^i^^ we extracted ki = 3.68, K2 — 0.66, a = 0.068 
which should be valid for a clean sample (single crystal). 
In these works the strength of the Josephson coupling, 
7, is not estimated. In the dirty limit, the compound is 
definitely a type-II material. Substituting former values 
into Eq. (|lip we obtain 5"^ /m « 0.59. On the other hand, 
from Refs.-i^'i^ we obtain Tc2 = 0.56Tci, and 5 = 1.33. 
For usual Mg^^B2 we have 7 — 0.4, m = 1 and a = 
0.844. For the Mgi"B2 we found 7 = 0.28. For irradiated 
MgB2 samples a mass ratio m k lA has been observed 
together with a = 0.059^^ Therefore, in the remainder 
of the paper we will not restrict ourselves just to the 
particular values of the parameters, but rather focus on 
new physics between the two types of superconductivity 
and its manifestations. 



A. Magnetic vs. Josephson coupling 

As mentioned above, one of the most fascinating prop- 
erties of a two-band system is the possible appearance 
of fractional and fractional states. We argued that those 
could be observed in experiment through their magnetic 



response. Using the preceding theoretical formalism, we 
can now calculate the magnetic field in and around the 
sample, in response to the applied magnetic field. In Fig. 
I18l we show the magnetic field profile emanating from the 
mesoscopic superconducting disk in the fractional (0,1) 
and (1,2) vortex state. The asymmetry induced by the 
Josephson coupling is clearly visible and can be directly 
imaged in magnetic force microscopy or scanning Hall 
probe magnctomctry experiments. Additionally, the in- 
tegrated magnetic field from such measurements will re- 
veal the fractional flux carried by these states of non- 
integer total vorticity. 

In Fig. [TH] we show the effect of both Josephson and 
magnetic coupling on the free energy and the vortex 
states. Fig. [T9la) shows the energy landscape when su- 
perconducting condensates are decoupled. As could be 
expected, all (L, L + 1) fractional states are stabilized 
in the ground state, at intermediate fields between inte- 
ger flux L and L + 1 states. (L, L + 2) states also exist, 
but have significantly higher energy. The same holds for 
(i + 1, L) fractional states, having far higher energy than 
the corresponding {L,L + 1) state. This can be entirely 
inverted for a different choice of parameters, particularly 
a, which determines the relative coherence lengths and 
consequently the ratio of the vortex energy in the two 
bands. In the present calculation, we therefore omit the 
curves corresponding to (L + 1, L) states. 

In Fig. [T9Fb) the Josephson coupling is added. This 
directly results in stabilization of the integer flux states 
in the ground state, and fractional ones have much higher 
energy. Generally, the (L, L + n) energy increases further 
and energy levels follow each other as magnetic field and 
n are increased. 

In Fig. Iigf c) we introduced the magnetic coupling, in 
absence of the Josephson one. This broadens the stabil- 
ity intervals of all vortex states - integer and fractional 
- as a consequence of the magnetic screening which low- 
ers the effectively experienced field by the sample and 
enhances superconductivity. However, the magnetic re- 
sponse of the sample is generally of oblate shape (due to 
the symmetry of the disk), and asymmetric states are less 
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FIG. 19: (Color online) Free energy as a function of magnetic 
flux for found integer and fractional states with L < 4. A 
linear background is subtracted from all curves to enhance 
readability (as indicated in the labels). The colors of the 
curves correspond to the vorticity, and also indicate the com- 
binations involved in the (two-colored) fractional states. The 
ground state is indicated by the color-coded shaded areas be- 
low the curves. Parameters of the sample are Rd = 7.5^0, 
T = 0.2, Tcr = 0.5, a = 0.3, and m = 1. 
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FIG. 20: (Color online) The flux quantization in a two gap 
superconductor bands only coupled by the magnetic field, i^r 
is the flux measured through the contour F, defined as the 
contour where the supercurrent equals zero, as illustrated by 
the white dashed line in the inset. The thick solid grey line 
in the graph depicts the ground-state as a function of applied 
field. 



favorable than in case (a) . For that reason, pronouncedly 
asymmetric (0,1) and (1,2) states are not present in the 
ground state. However, higher fractional states can be 
found in the ground state since the ring-like arrangement 
of vortices and their fields in both condensates enhance 
each other, and also follovir better the overall symme- 
try of the stray magnetic field. Nevertheless, Joscphson 
coupling is still able to completely remove the fractional 
states from the ground state, as shown in Fig. [TWd). 
Arguably, at higher vorticity fractional states could ap- 
pear in the ground state, since the order parameters in 
the two bands will become increasingly similar with in- 
creasing vorticity. 

In Fig. [5D]we illustrate the fiux-quantization in a two 
gap superconducting disk. As is already known, the fiux 
in single-band mcsoscopic samples is not quantized, but 
it always is within a contour determined by the zero cur- 
rent. Therefore we compute the fiux ^r penetrating the 
two-band sample through a contour F inside the sample 
on which j,, =0. The result is plotted in the presence of 
magnetic coupling between the bands to still have some 
fractional vortex states left in the ground state. By plot- 
ting the whole stability regions of the states, we noticed 
that the fractional flux decreases with increasing applied 
magnetic field (i.e. the applied fiux (/>). We pinpoint this 
effect in addition to the findings of Ref J^, where the au- 
thors found the decrease of flux through contour 7 with 
increasing temperature. The reason for the change of the 
fractional fiux with applied field or temperature is that 
one of the condensates always depletes faster than the 
other (in the present case, the second one). At sufficiently 
high field or temperature, only one band superconducts. 
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and the flux through contour T changes towards its quan- 
tized level in the surviving band. In the present case, the 
first band is stronger, and the fractional flux decreases 
towards vorticity in the first band, i.e. L = 1. If the con- 
sidered state was a (2, 1) one, the fractional flux would 
increase towards L = 2 level with increasing field or tem- 
perature. 



B. Magnetization curves 

We define the magnetization as: M = j^ J AHzdV, 
where AH represents the induced magnetic field. In Fig. 
HHwe plot the M{H)-\oop in a MgB2 disk with radius 
Rd = 10^0 and thickness d = ^ at T = 0, taking the pa- 
rameters of Ref ji^, i.e. rei = 3.68, K2 = 0.66, a = 0.068. 
Although the second band is conventionally type-I, the 
shape of the M {H)-loop suggests that the whole system 
still behaves like a type-II superconductor, i.e. there is 
no indication of a type 1.5 superconductivity reported by 
Moshchalkov et alJ^ We additionally plotted the magne- 
tization in the absence of Josephson coupling, which also 
does not show any qualitative deviation of type-II be- 
havior. The only influence of the Josephson coupling is 
an apparent increase of critical field and a stronger mag- 
netic response, which is a direct consequence of currents 
being strengthened by coupling. However, when we take 
lower values of the GL parameters, for example ki = 1 
and K2 = 0.2, we find, as shown in Fig. [33] a behavior of 
the magnetization versus field that is neither type-II nor 
type-I like. This state is characterized by a steep drop of 
the magnetization at a field close to the thermodynami- 

cal critical field He of the second condensate {He ) . Due 
to finite demagnetization effects, characteristic of type-I 

(2) 

samples, the transition is at a lower field than He ■ At 
this transition field, superconductivity ceases in the sec- 



ond gap, and the magnetization undergoes a steep drop. 
The origin of this effect is clearly visible in the figure, 
where also the mean Cooper pair density in the sample 
is plotted - the magnetization drop coincides with the de- 
pletion of the second band. Beyond the transition field, 
the flux continues to enter the sample gradually, exhibit- 
ing the type-II mixed state of the flrst condensate, and 
the overall behavior of magnetization can be treated as a 
superposition of type-I (steep drop) and type-II (gradual 
decrease) behavior of the two condensates, each being of 
different type. The influence of the Josephson coupling 
is also striking, as it smoothes out the drop in the mag- 
netization: the second condensate still depletes but at a 
slower rate due to the exchange of Cooper pairs with the 
flrst condensate. The slope of the decrease of the mean 
Cooper pair density in the second condensate still seems 
to match the slope of the drop in magnetization with 
a remarkable accuracy, although the transition becomes 
less abrupt and more reminiscent of a type-I intermediate 
state with bundles of flux penetrating the sample. 

In Fig. [33] and mi we demonstrate the influence of the 
ratio of the coherence lengths, a on the magnetic behav- 
ior of the sample. In Eq. p4)) . the supercurrents due to 
the second condensate have a prefactor of a/K^. There- 
fore, it is not K2 alone which determines the type of the 
band and an effective Kj = k-iI \foi can be introduced. 
Essentially, this suggests that the self induced fleld should 
be proportional to a, and the results in Fig. [331 are sup- 
portive of this. The level of magnetization to which the 
sample jumps after the magnetization drop should be in- 
dependent of a, since the second condensate is depleted 
there. However it turns out that the larger the magneti- 
zation is before the drop, the lower it becomes after the 
drop. This follows from the fact that, before the magne- 
tization drop, the second condensate is able to provide a 
better screening of the magnetic fleld for the flrst conden- 
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FIG. 21: (Color online) M{H) loops obtained by sweeping up 
and down the magnetic field for a MgB2 single-crystal disk of 
radius Rd = lO^io. 




FIG. 22: (Color online) Solid curves represent magnetization 
as a function of applied magnetic field for 7 = (black) and 
7 — 0.1 (red), and parameters of the sample ki = 1, K2 = 0.2, 
a = 0.1, Rd = 10^10 and T — 0. Dashed curves represent the 
mean Cooper pair density in the second condensate. 
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FIG. 23: (Color online) Magnetization of the sample versus 
the applied field for different values of the ratio between the 
coherence lengths in two condensates a. Other parameters 
of the sample are 7 = Ki = 1, reg = 0.2, Rn ~ lOCio and 

r = o. 
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FIG. 24: (Color online) Idem as Fig. [23l but for finite Joseph- 
son coupling 7 = 0.1. 



sate when a is larger, but when the second condensate 
ceases, the first condensate experiences a large difference 
in the felt magnetic field, which in turn allows for a larger 
flux penetration and thus a lower diamagnetic response 
of the sample. 

It should be stressed here that the field at which the 
second condensate depletes and the magnetization drops 
is also influenced by a. The smaller the coherence length 
of the second condensate (higher a), the smaller the tran- 
sition field. This difference is even more prominent in the 
presence of Josephson coupling. In other words, the ap- 
parent demagnetization effect in the type-I part of the 
magnetization curve is clearly influenced by the parame- 
ter a. 



Finally, in Fig. [25]we also show the magnetization cor- 
responding to the sweep-down of the applied magnetic 
field. At the point where the second condensate revives 
the magnetization jumps up, since there the type-I con- 
densate contributes to the diamagnetic signal. This jump 
is less abrupt when 7 is non-zero. 



C. Type-I. X vortex states 

It is already known that vortices repel each other in 
type-II superconductors, form Abrikosov lattice in bulk 
samples, but are compressed into geometry dependent 
multi-vortex states and even giant vortex states in meso- 
scopic superconductors. In type-I samples however, fiux 
penetrates the sample in the form of lamellae or tubular 
flux bundles»2^ As we showed above, two band supercon- 
ductors can show a bit of both behaviors, called type- 
1.5 superconductivity by Moshchalkov et al.M- While we 
demonstrated the type I.x behavior through the magne- 
tization loops, they mainly discussed the vortex-vortex 
interaction in two-band superconductors, claiming that 
it should be short range repulsive and long range attrac- 
tive. This of course assumes integer fiux vortex states, 
or strong Josephson coupling in our model. However, as 
we have seen above, a plethora of other vortex states are 
possible, not all with integer flux. Therefore, the vortex- 
vortex interaction should be discussed separately within 
bands (intra-band) , and separately between them (inter- 
band). While leaving the detailed analysis for the future, 
we here show several prime examples of vortex states 
that can be found in two-band mesoscopic disks, that 
show type-I. X behavior (however different from RefJ^). 
For example, we take the disk with radius lO^io with 
parameters of the condensates ki = 1, K2 = 0.2 and 
a = 0.3. As shown in Fig. I^TTa-b). at larger applied 
field providing 30 flux quanta through the sample, we 
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FIG. 25: (Color online) Full sweep up and down M{H) loops 
with and without coupling. Taken parameters are a — 0.3 
Ki = 1, K2 = 0.2, Rd = 10^10 and T = 0. 
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FIG. 26: (Color online) M{H) for all stable vortex states in 
a mesoscopic type-I.x superconducting disk. 



found a L = 24 vortex state. Due to absence of Joseph- 
son coupling, two allotropic modifications of the vortex 
state were found possible. In both, the vortices in type-I 
condensate merge into a single domain, either a ring do- 
main (a), or a giant-bubble (b), which is typical for type-I 
samples. Actually they are still single vortices, but have 
huge overlap, mimicking a normal domain. Vortices in 
the type-II condensate remain separate however, but are 
forced by magnetic coupling to obey the symmetry of the 
intermediate state of the type-I condensate. They there- 
fore form more or less a conventional multi-shell state 
under the giant vortex in (b) , but are forced to make an 
unconventional state with dense shells within the ring do- 
main in (a) . This type-TII competition is even more pro- 
nounced in the presence of Josephson coupling, as shown 
in Fig. [2Wc). In the first condensate, one can see three 
giant L = 2 vortices and two L = 3 multi-vortex clus- 
ters. This is also the situation in the second condensate, 
however the multi-vortices are now even closer and also 
overlap more, mimicking perfect giant vortices. Due to 
the Josephson coupling, both condensates influence each 
other; as a consequence, the type-I intermediate state is 
forced to split into as many bubbles as possible, and vor- 
tices in type-II condensate must group into those bubbles. 
As a result, a multi-vortex of multi-vortices is formed, 
clearly a signature of type-I.x behavior. 



IV. CONCLUSIONS 

In summary, we presented a theoretical Ginzburg- 
Landau (GL) study of the superconducting state of 
two-band mesoscopic disks, where both the influence of 
Josephson and of magnetic coupling between the super- 
conducting bands are discussed. In cases when screening 
of the magnetic field can be neglected, we found the de- 
pendence of the size of the vortex core on the strength of 
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FIG. 27: (Color online) Cooper pair density contour plots 
of the first (left) and second (right) condensate. ki — 1, 
K2 = 0.2, a = 0.3, Rd = 15^10- In (a) and (b) (ji/fjio = 30, 
L = 24, 7 = 0. In (c) 0/0o = 15, L = 12, 7 = 0.02. 



the Josephson coupling and showed that it generally has 
an influence opposite to the one of temperature. In limit- 
ing cases, our numerical flndings agree well with analytic 
expressions available in literature. We also found a flt- 
ting function, which gives an excellent estimate of the size 
of the vortex core as a function of the size of the meso- 
scopic disk. In our further analysis of the vortex states, 
we focussed mainly on exotic, fractional states, where 
two condensates comprise different number of vortices 
and the apparent total vorticity of the sample is frac- 
tional. We reported asymmetric vortex states following 
from competing interactions in the two-band mesoscopic 
system, and showed how some states can be manipulated 
by e.g. coupling between the bands. We indicate how 
such states can be experimentally observed. Fractional 
states can even be found in the ground state, but typi- 
cally far from the S/N boundary. We give the expression 
for the upper critical field of a two-band mesoscopic disk 
as a function of temperature, which is similar to analytic 
estimations for bulk, however scaled to its value at zero 
temperature for zero coupling between the condensates. 

When magnetic screening and coupling between the 
bands is included in the simulations, we characterized 
the response of the sample through the competition of 
the GL parameters of the two-bands (with special atten- 
tion to the case when one band is typc-II and the other 
is typc-I). However, we show that this is insufficient, and 
that Josephson coupling and the ratio of the coherence 
lengths in the two bands also play an important role. 
Although we did not find evidence for type-1.5 supercon- 
ductivity in clean MgB2 disks, we did find its manifes- 
tation for a different choice of relevant parameters. The 
magnetization vs. applied field shows a distinct jump at 
the field where type-I condensate ceases, and the overall 
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shape of the curve can surely be characterized as type- 
I.x Uke. This is also evident in the found vortex states in 
the latter case, which arc a combination of single vortices 
and lamellar domains. 
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